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Abstract 

Metabolomics is providing new dimensions into understanding the intracellular adaptive responses in plants to external 
stimuli. In this study, a multi-technology-metabolomic approach was used to investigate the effect of the fungal sterol, 
ergosterol, on the metabolome of cultured tobacco cells. Cell suspensions were treated with different concentrations (0- 
1000 nM) of ergosterol and incubated for different time periods (0-24 h). Intracellular metabolites were extracted with two 
methods: a selective dispersive liquid-liquid micro-extraction and a general methanol extraction. Chromatographic 
techniques (GC-FID, GC-MS, GCxGC-TOF-MS, UHPLC-MS) and 'H NMR spectroscopy were used for quantitative and 
qualitative analyses. Multivariate data analyses (PCA and OPLS-DA models) were used to extract interpretable information 
from the multidimensional data generated from the analytical techniques. The results showed that ergosterol triggered 
differential changes in the metabolome of the cells, leading to variation in the biosynthesis of secondary metabolites. PCA 
scores plots revealed dose- and time-dependent metabolic variations, with optimal treatment conditions being found to be 
300 nM ergosterol and an 18 h incubation period. The observed ergosterol-induced metabolic changes were correlated 
with changes in defence-related metabolites. The 'defensome' involved increases in terpenoid metabolites with five 
antimicrobial compounds (the bicyclic sesquiterpenoid phytoalexins: phytuberin, solavetivone, capsidiol, lubimin and 
rishitin) and other metabolites (abscisic acid and phytosterols) putatively identified. In addition, various phenylpropanoid 
precursors, cinnamic acid derivatives and - conjugates, coumarins and lignin monomers were annotated. These annotated 
metabolites revealed a dynamic reprogramming of metabolic networks that are functionally correlated, with a high 
complexity in their regulation. 
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Introduction 

Metabolomics is a holistic qualitative and quantitative analysis 
of all metabolites present within a biological system under specific 
conditions [1-4]. Metabolomics differs from the classical or 
traditional targeted phytochemical analysis in various fundamental 
aspects such as being a data-driven approach with predictive 
power that aims to assess all measurable metabolites without any 
pre-conception or pre-selection. In order to attain this goal, 
advanced analytical tools that provide high degrees of sensitivity, 
selectivity and reproducibility are required [3,5-7] . Metabolomics 
is viewed as a complementary technique to other functional '-omics' 
approaches such as transcriptomics and proteomics. The integra- 
tion of the these technologies contributes to a systems biology 
overview [8] , providing a holistic understanding of the organisa- 
tion principle of cellular functions at different levels, and ways of 
monitoring all biological processes operating as an integrated 
system [4,9-1 1]. Moreover, metabolomics as a post-genomics tool 
is often regarded as offering distinct advantages when compared to 
other 'omics' technologies. This point of view is based on the fact 
that changes in the transcriptome or proteome do not always 
correlate to biochemical phenotypes [4,7,12,13]. 



The holistic analysis of the metabolome, with its complex/ 
divergent physico-chemical properties and dynamic molecular 
composition, requires a wide range of chemistries. Hence, all 
metabolomic analyses are like a snapshot (or point-in-time- 
chemistry) of a biological system (cell, tissue or whole organism), 
showing which metabolites are present and the levels at a given 
time point and under specific physiological conditions 
[1,3,7,11,14]. Different strategies and a range of analytical 
techniques have thus been developed for different metabolomic 
analyses; and the usage of parallel analytical platforms can provide 
a wide coverage of the metabolome under study, additional 
information or confirmation for a putatively identified metabolite 
[3,14-16]. 

In plant research, metabolomic approaches are increasingly 
being used for various studies including linking genotype and 
biochemical phenotype, silent phenotypic mutations, metabolic 
pathway studies, and abiotic- and biotic stresses, including plant : 
pathogen interactions [2,4,9,10,13,17-20]. Plants are continuously 
threatened by a wide range of pathogens or abiotic stresses and, 
for protection, all plants possess well-established resistance 
mechanisms developed through evolution. The innate immune 
system of plants can be divided into two fundamental components: 
protection and defence [21,22]. Protection is a static (passive) 
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phenomenon and involves both structural barriers and pre-formed 
inhibitors. These protective mechanisms prevent or attenuate 
invasion by potential attackers [21,23-25]. On the other hand, 
defence is an inducible, dynamic (active) phenomenon and occurs 
only when the host and the pathogen have made physiological 
contact [21,22]. The final outcome of the plant : pathogen 
interaction depends ultimately on the balance between the ability 
of the pathogen to suppress and overcome the plant's immune 
responses, and the capacity of the plant to recognise the pathogen 
and activate effective defences [26,27]. Recognition involves the 
binding of microbe-associated molecular pattern molecules 
(MAMPs) to pattern recognition receptors (PRRs) in the host. A 
detailed understanding of plant resistance mechanisms (protective 
and inducible response) opens up possibilities of developing 
strategies to increase stress/disease resistance in plants, noting 
that the latter are an essential part of human life (foods, medicines, 
industrial raw materials, etc.) and a vital component of the whole 
ecosystem. 

Among the phytopathogens are fungi that continue to pose a 
significant threat to crop production and subsequently to food 
supply. The early stages of plant : pathogenic fungal interactions 
are mostly mediated by fungal MAMP molecules such as the 
evolutionary conserved ergosterol; a 5,7-diene sterol. Ergosterol 
does not occur naturally in plants and differs structurally from the 
phytosterols, having two additional double bonds (at positions C7- 
C8 and C22-C23) and a methyl group at C24 of the side chain. It 
is thus recognised by such cells as 'non-self [28], thereby acting as 
a MAMP molecule in the lipophilic class of biotic elicitors [29]. 
The effect of ergosterol on plant secondary metabolism has not 
been thoroughly investigated [30-32]. Following our previous 
report [33], here we present a multi-platform metabolomics-based 
elucidation and analysis of changes in the metabolism of tobacco 
(Mcotiana tabacum) cells following ergosterol treatment. The results 
significandy contribute to a more comprehensive understanding of 
the defence-related phenotypic status of plants. 

Materials and Methods 

Chemicals and Reagents 

All chemicals used were of analytical and ultra-pure LC-MS 
grade quality, and organic solvents mostly included chloroform 
(Labscan, Poland), methanol (Labscan, Poland) and acetone 
(Associated Chemical Enterprises, South Africa). All equipment 
was sterilised prior to use and cell treatment was carried out under 
sterile conditions. 

Cell Culture and Elicitation of Tobacco Cell Suspensions 

Cell suspensions are used in metabolomics studies as systems 
that allow reproducibility and controllability of conditions [34— 
36], and Mcotiana tabacum cv. Samsun cell suspensions, cultivated as 
previously described [37,38] were used. For elicitation of cells, a 
stock solution of 4.54 mM ergosterol, C28H44O, (Sigma, USA), 
was prepared in acetone. The final concentration of acetone in the 
treated suspensions did not exceed 0.05% and had no discernable 
effect on the observed responses. 

Three days after subculture, cells were treated by adding 
specific volumes of the stock solution of ergosterol to aliquots of 
cells suspensions with continuous rotation at 80 rpm and 25°C. 
For concentration studies, ergosterol was added to final concen- 
trations of 0-1000 nM, with an 18 h incubation period while non- 
treated cell suspensions were used as negative controls. For the 
time study, cell suspensions were treated with 300 nM ergosterol 
and incubated for 0-24 h with a non-treated sample, incubated for 
24 h, included as a control. One of the key points in metabolomic 



studies lies in differentiating true biological variation from 
technical variation. Three independent biological repeats were 
therefore conducted and analysed. 

Cell Viability Determination 

In order to determine if possible secondary responses due to 
treatment-induced cell death occurred, a cell viability assay was 
performed based on the reduction of 3,5-triphenyltetrazolium 
chloride (TTC) to red triphenylformazan (TPF) as previously 
described [33,39]. 

Metabolite Extraction 

Extraction methods for metabolomic studies depend on the 
chemical and physical properties of the target metabolites, the 
biochemical composition of the system from which the analytes are 
to be extracted and the properties of the solvents to be used 
[15,40]. The main extraction method developed to extract these 
metabolites was the dispersive liquid-liquid microextraction 
(DLLME) technique [41] using methanol and chloroform as 
disperser and extraction solvents respectively [33]. Briefly, 
DLLME was carried out as follows: after induction with ergosterol 
and incubation, the excess media was filtered using a Buchner 
funnel and the cells collected. Two grams of cells were re- 
suspended in 20 mL of methanol and homogenised using an Ultra 
Turrax homogenizer. The homogenates were centrifuged at 
5525 xg for 7 min at 25°C. The supernatants were placed into 
clean 50 mL round-bottom flasks, and the methanol evaporated to 
1 mL at 50°C using a Buchi Rotavapor. One mL aliquots of the 
resulting crude aqueous extracts were transferred to glass 
centrifuge tubes and 200 |0,L of chloroform and 100 |JL of 
methanol were added. The mixture was vortexed for 30 s and 
centrifuged at 9 000 xg for 6 min at 25°C. Using a syringe, the 
bottom chloroform extraction layer was removed and filtered 
through 0.22 |J,m filters (Millipore, USA), and placed in glass vials. 
The chloroform extracts were kept at — 20°C until further analysis. 

A general methanol extraction was also used with the aim of 
extracting more polar metabolites. The first steps of this methanol 
extraction, from the induction to the Buchi Rotavapor methanol- 
evaporation step, were carried out as described above in the 
DLLME method. The 1 mL crude aqueous extract (following the 
evaporation of methanol) was placed into clean 2 mL microcen- 
trifuge tubes and evaporated to dryness under vacuum using a 
SpeedVac centrifugal evaporator at 55°C. To the dry extract, 
100 |XL of methanol and 100 jlL MilliQ water were added. The 
mixture was vortexed until all content was dissolved. The 50% 
methanol extracts were then filtered through 0.22 (Xm filters and 
placed into glass vials for chromatographic analyses. The 
methanol extracts were kept at — 20°C until further analysis. 

Following the extraction of the metabolites by the two methods, 
the extracts were further analysed, qualitatively and quantitatively, 
using various chromatographic techniques, mass spectrometry, 
and one-dimensional proton-nuclear magnetic resonance spec- 
troscopy. Hence, to ensure that the data acquired from these 
analytical techniques were acceptable for further multivariate data 
analysis, three technical repeats for each biological repeat were 
analysed, for all experiments. 

NMR Spectroscopy Analyses of Chloroform Extracts 

One-dimensional proton nuclear magnetic resonance spectros- 
copy (ID 'H NMR) analyses were performed using "1 PULSE" 
Fourier-transform (FT)-NMR experiments. All spectra were 
measured using a Bruker AVIII-400 NMR spectrometer, operat- 
ing at a proton NMR frequency of 400.17 MHz at 296.8 K. For 
each sample, 64 scans (number of scans, NS) were recorded using 
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a 9.30 \ls (60°) pulse (PI, pulse width/pw), 8.22 kHz spectral 
width (SWH), and 1.00 s relaxation delay (Dl/dl), and 0.125 Hz/ 
point (FIDRES/res). The acquisition time (AQ/at) was set at 
3.98 s and the free-induction decays (FIDs) were Fourier 
transformed with a line broadening (LB) factor of 0.30 Hz. For 
quantitative analysis, peak height was used. The CDCl :i solvent 
used contained the internal standard, trimethylsilyl propionate 
(TSP). The spectra were referenced to both residual solvent signal 
of CHC1 3 (7.240 ppm) and TSP at 0.000 ppm. 

Gas Liquid Chromatography Coupled to Flame lonisation 
Detection (GC-FID) and Electron-Impact-Mass 
Spectrometry (GC-EI-MS) Analyses 

DLLME/ chloroform extracts were analysed on a Shimadzu- 17- 
10A GC instrument (Shimadzu, Kyoto, Japan), equipped with a 
30 mx().25 mmx0.25 (Xm Zebron ZB-1MS column (Phenom- 
enex, California, USA). A 2 |lL aliquot was injected into the GC- 
FID instrument in splitless mode. The injection and detector 
temperatures were set to 250°C and 350°C respectively. The 
overall programmed-temperature GC run included zones with 
different temperature gradient (multi-linear programming) with 
initial and final isothermal zones. The initial GC oven temperature 
was 80°C for 4 min, increased by 25°C min -1 to 240°C, then 
increased to 300°C at a rate of 5°C min - 1 and held at 300°C for 
3.60 min. The total run time was 26.00 min. The flow rate was set 
to 1 mL min - , with nitrogen as the carrier gas. The detection 
mode in this GC approach was flame ionisation detection (FID). 
The solvent (chloroform) blank was injected after baking the 
column and cleaning the injection port on the GC instrument to 
determine if any peaks were attributable to the solvent (no 
significant peak was present apart from the solvent peak). This was 
done prior to the analysis of each biological replicate. 

GC-MS was generally used for identification purposes as MS 
detection gives structural information about the analytes. In this 
experiment, a Shimadzu GC-MS-QP2010 (Shimadzu, Kyoto, 
Japan) with an electron impact (EI) ion source and a quadrupole 
mass analyser was employed. The column used was a ZB1-MS 
non-polar column (Phenomenex, California, USA). GC condi- 
tions/ parameters were set as described above. The injection mode 
was set to a splitless; and 1 |0,L of the chloroform sample was 
manually injected. The MS was operated using a positive EI at 
70 eV. Helium was used as the carrier gas, and the column flow 
rate set at 1 mL min . The ion source temperature was set at 
200°C, and the detector voltage was set at 0.5 kV to prevent the 
oversaturation of the detector since the chloroform extracts were 
concentrated. The solvent cut time for the MS was set at 5 min 
and EI scan mode was used for identification covering the range of 
50-500 m/z. 

Comprehensive 2D GC Coupled to TOF-Mass 
Spectrometry Analyses 

In the GC xGC technique the two columns are coupled by a 
special interface (modulator) that is capable of sampling the 
effluent from the first column (first dimension) and periodically 
introducing it to the second column (second dimension) in a 
manner that preserves the original first dimension separation [42- 
44]. The DLLME chloroform extracts were analysed using a 
Pegasus 4D GCxGC-TOF-MS instrument from LECO Corpo- 
ration (St. Joseph, MI, USA). Samples were run in two- 
dimensional mode. The sample injection volume was 1 |jL with 
helium as the carrier gas at a flow rate of 1.00 mL min The 
inlet temperature was set at 250°C and a splitless mode injection 
type was used. The first dimension column was a mid-polar Restek 



Rxi-17 S matrix column (15 m, 360 |J,m i.d., 0.25 |Jni d.f.) and the 
temperature profile was 10.00°C min -1 from 50 to 325.00°C with 
a hold for 1.00 min at 50°C and for 10.00 min at 325.00°C. The 
second dimension column was a non-polar Restek Rtx-5 matrix 
column (0.990 m, 180 |0.m i.d., 0.20 |Im d.f.) and the temperature 
profile was 10.00°C min 1 from 65 to 340°C with a hold for 
1.00 min at 65°C and for 3.00 min at 340°C. The first and second 
dimensional columns were linked by a modulator, and the 
modulation timing was set at 5.00 s. The detector was a time-of- 
flight (TOF)-mass spectrometer, and electron impact (EI) was used 
as the ionisation mode. The ion source temperature was set at 
200°C while the detector voltage was set at 1750 V and the 
electron energy at —70 V. The mass range collected was 40- 
550 Da with an acquisition rate of 100 spectra s 

Ultra High Performance Liquid Chromatography-High 
Definition Mass Spectrometry (UHPLC-HDMS) Analyses 

UHPLC and High Definition Mass Spectrometry (UHPLC- 
HDMS) analyses were performed on a Waters Acquity UHPLC 
coupled in tandem to a Waters photodiode array (PDA) detector 
and a SYNAPT Gl HDMS QTOF mass spectrometer (Waters, 
Manchester, UK). The combination of UHPLC with a PDA and 
QTOF detectors allow for the tentative identification of 
compound classes, and the possibility to obtain a more conclusive 
identification of separated compounds utilising mass spectrometric 
measurements and software tools. The HDMS technique (Waters, 
Manchester, UK) provides much structural information, and 
particularly high mass measurement accuracy, typically better 
than 5 mDa, which is a key factor in compound identification in 
metabolomics. This HDMS technology expands beyond the 
conventional MS due to its capabilities such as extra-dimensional 
high efficiency ion mobility, improved analytical peak capacity and 
selectivity, enhanced sensitivity, efficient sampling of ions, and 
elimination of neutral contaminants [45-47]. 

Mass accuracy was obtained by calibration of the instrument 
over the mass range of the analytical method and was further 
enhanced by the use of a reference mass channel generated by the 
use of a lock spray interface and leucine enkephalin as calibrant. 
Any fluctuations in the accuracy of mass measurements, caused by 
temperature fluctuations and electronic noise, were automatically 
corrected if needed. This design allowed for extended run times as 
well as batch analysis spreading over several days without 
compromising mass accuracy. 

Chromatographic separation of both the chloroform and 
methanol extracts was done utilising a Waters CSH C18 Acquity 
column (150 mmx2.1 mm, 1.7 |0,m) thermostatted at 60°C. A 
binary solvent mixture was used consisting of water (eluent A) 
containing 10 mM formic acid (pH of 2.3) and acetonitrile (Romil 
pure chemistry, UK) (eluent B). The initial conditions were 95% A 
at a flow rate of 0.4 mL min - 1 and were kept constant for 2 min. 
A gradient was introduced to change the chromatographic 
conditions to 5% A at 22 min. The conditions were kept constant 
for 3 min to flush the analytical column whereafter the column 
was returned to initial conditions at 27 min and allowed to 
equilibrate for 3 min. The run time was 30 min and the injection 
volume was 5 |ll. Each sample was analysed in triplicate, to 
account for any analytical variability. The PDA detector was 
scanned between 200 and 500 nm (1.2 nm resolution) and 
collecting 20 spectra s 

The SYNAPT Gl mass spectrometer was used in V-optics and 
operated in electrospray ionisation (ESI) mode to detect the 
compounds of interest. Leucine enkephalin (50 pg mL ') was used 
as reference calibrant to obtain typical mass accuracies between 1 
and 3 mDa. The mass spectrometer was operated in both positive 
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and negative mode with a capillary voltage of 2.5 kV, the sampling 
cone at 17 V and the extraction cone at 4 V. The scan time was 
0.1 s covering the 100-1000 Da mass range. The source 
temperature was 120°C and the desolvation temperature was set 
at 450°C. Nitrogen gas was used as the nebulisation gas at a flow 
rate of 800 L h _1 . The software used to control the hyphenated 
system and perform all data manipulation was MassLynx™ 4.1 
(Waters Corporation, USA). 

Multivariate Data Analysis 

In comparison to traditional univariate statistical methods, the 
MVDA models are well suited to provide ways of handling 
confounding and covariance patterns (both within and between 
variables), which are found in complex and multi-dimensional data 
sets from metabolomics studies [48,49]. Principal component 
analysis (PGA), an unsupervised multivariate linear model, and 
orthogonal projection to latent structures-discriminant analysis 
(OPLS-DA), a supervised model, were used for data analysis. Only 
GC-FID, UHPLC-MS and NMR data were MVDA-modelled, 
generating sufficient and conclusive information which was 
supplemented by visual inspection of chromatograms and NMR 
spectra. 

The quality of the models was evaluated based on some model 
diagnostic tools largely used in metabolomic studies, namely the 
cumulative modelled variation in matrix X, irX(cum) or the 
goodness-of-fit parameter, the proportion of the variance of the 
response variable that is explained by the model, .ff^cum), and 
the fraction of the total variation of matrix X that can be predicted 
by the extracted components, Qf(cum) known also as predictive 
ability parameter. The values of these diagnostic parameters must 
be close to 1.0 for a robust mathematical model with a reliable 
predictive accuracy [50,51]. 

NMR Data Analysis 

'H NMR spectra were automatically reduced to ASCII files 
using AMIX software (version 3.7, Bruker Biospin, Germany). 
Spectral intensities were scaled to total intensity and reduced to 
integrated regions of equal width (0.04 ppm) corresponding to the 
region of 8 0.1-10.00. The region of 8 7.3 was excluded from the 
analysis because of the residual signal of chloroform. The data 
matrix obtained was exported to the SIMCA-P software version 
13.0 (Umetrics, Umea, Sweden) for PCA modelling using the 
Pareto scaling method. The Panto scaling involves dividing the 
variables by the square root of their standard deviations [49,52- 
54] . PCA scores and loadings plots were used to explain variations 
in the samples. 

GC-FID Data Analysis 

For GC-FID data analysis, a data matrix was generated from 
the retention time (Rt) and peak area chromatographic data (from 
the Shimadzu-17-10A GC instrument). This dataset was exported 
into the SIMCA-P software version 13.0 (Umetrics, Umea, 
Sweden) for PCA modelling. The PCA scores plot was used to 
depict a visual image of sample variations from a global view. The 
data were all iWfo-scaled so as to reduce the impact of noise and 
artefacts in the models, improving subsequently the models' 
predictive ability. 

UHPLC-MS Data Analysis 

For UHPLC-ESI-MS data, both PCA and OPLS-DA models 
were used. ESI positive and negative raw data were extracted 
using MassLynx™ XS software and analysed with Marker- 
Lynx™ XS software (Waters Corporation, Mildford USA). The 



MarkerLynx software extracts the raw LC-MS data and produces 
a matrix of Rt-mlz variable pairs, with the ml z peak intensity for 
each sample. MarkerLynx software parameters were set to analyse 
the 2-26 min Rt range of the chromatogram, mass range 100— 
700 Da, mass tolerance 0.01 Da, mass window 0.05 Da and a Rt 
window of 0.20 min. The data matrix obtained from MarkerLynx 
processing was also exported to the SIMCA-P software for PCA 
and OPLS-DA modelling. In both models, the data were Pareto- 
scaled. PCA scores and loadings plots were used to explain 
variations in the samples. The OPLS-DA S-plot aided in 
biomarker/compound identification. The S-plot helps in explain- 
ing and visualising the covariance and correlation between ions 
and the modelled classes, allowing thus the mining of the data to 
extract metabolites that are statistically interesting compounds 
with potential biochemical significance [54]. Furthermore, the 
significance of the selected mass ions (variables) in OPLS-DA 
models was also assessed using the variable importance in 
projection, (VIP) plot (Umetrics, Umea, Sweden). 

Metabolite Annotation 

GC-MS. For GC-EI-MS the respective mass spectra of the 
extracted ion peaks identified by comparison of the total ion 
current (TIC) mass chromatograms were searched against the 
NIST (National Institute of Standards and Technology) and Wiley 
mass spectra libraries (2007 versions) for similarity matches. This 
spectral comparison provided putative empirical formulae and 
structures, which were further searched in databases such as the 
Dictionary of Natural Products (DNP) (www.dnp.chemnetbase. 
com) and ChemSpider (www.chemspider.com). For GCxGC- 
TOF-MS analyses, the Chroma-TOF software (LECO Corpora- 
tion, St. Joseph, MI, USA) was used for automated data analysis 
from preprocessing to the compound identification step. 

UHPLC-MS. Selected mass ions from S- and VIP plots were 
annotated using the Taverna workbench for PUTMEDID_LCMS 
Metabolite ID Workflows (www.taverna.org.uk). The Taverna 
workflows allow for integrated, automated and high-throughput 
annotation and putative metabolite identification from ESI-LC- 
MS metabolic data. The workflows consist of correlation analysis, 
metabolic feature annotation and metabolite annotation. Marker- 
Lynx XS data were exported into the Taverna workflows to be 
processed for metabolite annotation. The data matrix from 
MarkerLynx processing was firstly formatted to match the 
Taverna workbench requirements. Three main workflows formed 
the Taverna metabolite annotation procedure. In the workflow 1 
(List_CorrData), all the parameters were set at default, ion mode 
selected (as negative or positive depending on the ion mode in 
which the data where acquired) and Pearson correlation 
calculation was chosen as the correlation method. In the workflow 
2 (annotate_Massmatch), mass tolerance was set to 5 ppm, retention 
time range set as 120-1560 s (2-26 min) and all other parameters 
were set at default. These workflows allowed for grouping together 
ion peaks with similar features such as Rt, and annotating features 
with the type of ml z ion (molecular ion, isotope, adduct, others) 
believed to originate from the same compound. The elemental 
composition/molecular formula (MF) of each ml z ion is then 
automatically calculated. In workflow 3 (matchA'IF-MF) the 
calculated MF (from the workflow 2 output file) was then 
automatically compared and matched to the MF from a pre- 
defined reference file of metabolites (a customised chemical library 
based mainly on PlantCYC and AraCYC libraries). The selected 
MF were also manually searched against freely online databases 
such as DNP, Chemspider and METLIN libraries for more 
confidence in the annotated metabolites (Mi-level 2 annotation). 
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Results and Discussion 

The triphenyltetrazolium chloride-based cell viability assay 
[39], was performed to determine if treatment-induced cell death 
occurred under the experimental conditions. The results (not 
shown) indicated that treatment of the cell suspensions with 
ergosterol (0-1000 nM) did not lead to a significant cell death and 
that >91% of cells remained viable as estimated at 1000 nM, 
incubated for 18 h. The observed responses are thus due to the 
ergosterol treatment alone and possible secondary responses due to 
cell death can be excluded. 

In order to investigate the ergosterol-induced changes in the 
intracellular metabolome (as present in both semi-selective 
DLLME chloroform extracts and a general methanol extraction 
(ME) of tobacco cells, concentration- and time studies were 
conducted. These investigations served also to establish the 
optimal conditions for treatment of cells, which were found to 
be 300 nM ergosterol and an 1 8 h incubation period as discussed 
in this section. One-dimensional proton NMR spectroscopy and 
different chromatographic techniques (GC-FID, GCxGC-TOF- 
MS, UHPLC-MS) were used for the analyses. The combination of 
these analytical platforms provided a broader and more conclusive 
observation of the effect of ergosterol treatment on the metabolism 
of tobacco cells. 

One-Dimensional Nuclear Magnetic Resonance (1D 
NMR) Spectroscopic Analyses 

Chloroform extracts, obtained through the DLLME procedure, 
were analysed by ID H NMR to investigate metabolite variations 
due to different ergosterol concentration treatments. The NMR 
technology has been widely used as a fingerprinting tool in 
metabolomic analyses, with pattern recognition techniques such as 
PCA modelling [35,55,56]. NMR is a non-biased, non-destructive 
technique and one of the few analytical platforms that are rapid, 
reproducible, stable and require very littie or no sample 
preparation [1,5,57]. A 'H NMR spectrum of a plant extract is 
the result of the superposition of the 'H NMR spectra of all NMR- 
visible single compounds present therein. Hence, a proton NMR 
analysis would give a global view of all the metabolites (primary 
and secondary) in a sample, provided there are 'H NMR 
detectable [5,57]. 

Here, the DLLME samples were analysed using ID 'H NMR 
spectroscopy since the aim was to generate a global view or a 
fingerprint of the effect of different ergosterol concentrations on 
tobacco cells from a 'H NMR perspective. Figure 1 represents the 
results obtained from the ID 'H NMR analyses of the DLLME 
concentration study samples. 

The visual inspection of the overlaid NMR spectra of ergosterol- 
treated and non-treated control samples (Figure 1A) show 
treatment-related differences in various spectral regions: 80-83, 
83-85.8 and 85.8-810. These spectral regions are often known to 
be associated with amino acids, sugars and aromatic compounds, 
respectively [57]. However, other compounds also resonate in 
these regions. For instance, fatty acids and terpenes resonate in the 
0-3 ppm region and the sesquiterpenoids, rishtin and capidiol, 
also resonate in the 3-5.8 ppm spectral region [58]. 

The obtained NMR data were analysed by unsupervised 
MVDA modelling, PCA, to differentiate between ergosterol 
treatments, and to provide interpretable visualisation of the 
NMR spectra. PCA reduces the dimensionality of the data 
without much loss of information, and expresses the data in such a 
way as to identify and highlight the similarities and differences in 
systematic patterns and features of the data set. The PCA scores 
plot offers a visual image of sample variations from a global view 



and, being a non-parametric analysis, the generated model is 
independent of the user, hence unsupervised [48,49,59]. A five- 
component model, explaining 99.6% of the variance (with the 
accuracy of prediction of 0.991), was calculated. Scores plots PCI 
vs. PC 2 [i^X(cum) of 0.98, Q?(cum) of 0.969 and 95% confidence] 
and PCI vs. PC 3 [i^X(cum) of 0.91, Qf(cum) of 0.947 and 95% 
confidence] were constructed (Figures 1-B and -C), showing a 
differential clustering of the samples. The non-treated samples 
appear separate from the ergosterol-treated samples. The samples 
treated with 300 nM and 1000 nM ergosterol grouped together, 
with less variation between these two clusters. 

This NMR-based global view, as depicted by the PCA scores 
plots (Figure 1-B and -C), demonstrates that ergosterol treatment 
led to differential changes in the metabolite composition of the 
tobacco cells. The loadings plots (Figure 2) single out the 
discriminating variables responsible for sample clustering to be 
putatively compounds resonating in the 3-5.8 ppm spectral region 
(83.5 and 84.9 chemical shifts), and compounds abundantly in the 
amino acid/fatty acid spectral region (80-83 chemical shifts). 
Considering the low polarity of chloroform as solvent, it could be 
possible that the chemical shifts that contributed to the sample 
clustering are from sesquiterpenoids (in the 3-5.8 ppm spectral 
region) [58], fatty acids (in 0—3 ppm spectral region) and related 
precursor molecules [60-63]. However, no further definitive 
chemical elucidation of these NMR-detected compounds was 
carried out, which might require 2D NMR analyses to resolve any 
overlapping signals, or hyphenating NMR and LC to introduce a 
separation dimension [5,57]. Here, however, the NMR global 
view was sufficient to demonstrate ergosterol-induced metabolic 
changes. Furthermore, the inherent low sensitivity of NMR, its 
sensitivity to the chemical environment (pH, ionic strength, 
temperature, etc.) of the sample and the differential sensitivity of 
metabolites to the chemical environment hamper the quality of 
NMR analyses of complex samples [57]. 

Gas Chromatography-Flame lonisation Detector (GC-FID) 
Analyses 

Chloroform extracts, obtained through the DLLME procedure 
were analysed by GC-FID to investigate metabolite variations due 
to different ergosterol treatment conditions with regard to 
concentration and elicitation time. The samples were analysed in 
multi-linear programmed-temperature conditions to obtain a 
greater chromatographic separation. The resulting chromato- 
grams (Figures 3-A and -B) showed concentration- and time- 
dependent variations in terms of peak intensities and presence/ 
absence of peaks. The raw data were then preprocessed (peak 
alignment) and pre-treated (defining selected peaks) for further 
data analyses. 

The GC-FID chromatograms (Figures 3-A and -B) indicate that 
various compounds in the DLLME samples were separated and 
detected, showing treatment-related variations of metabolites with 
low polarity and molecular weight [64] , thus indicating changes in 
metabolite levels and expression, i.e. an altered metabolome. 

PCA was carried out to differentiate between ergosterol 
treatments. For the concentration study, a four-component model 
was generated and explained 94.7% of the variance [iTX(cum) of 
0.947 and Qf(cum) of 0.810], with the first two principal 
components (PCI and PC2) explaining 87.8%. A scores plot was 
constructed using PCI and PC 2 [R z X(cum) of 0.878 and Qf(cum) 
of 0.804; with 95% confidence level], showing cellular samples 
differentially clustered into five groups (Figure 4A). The clusters of 
50-, 150- and 300 nM ergosterol-treated samples are found to be 
grouped together, which implies a similarity in the GC-FID 
detected metabolite composition of these samples. The non- 
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Figure 1. Proton NMR spectra of extracts from ergosterol-treated cells. (A) Overlaid NMR spectra of DLLME samples representative of 
control, 300 nM and 1000 nM ergosterol-treated tobacco cells showing differences in proton shifts/signals due to ergosterol-induced changes across 
the 1-8 ppm region of the spectrum. The 'H NMR spectrum of ergosterol is at the bottom. Deuterated chloroform was used as the solvent; however, 
since deuteration is not 100%, the residual protons from chloroform (CHCI 3 ) give a singlet signal at 7.26 ppm. (B) and (C) PCA scores plots of the 'H 
NMR data from the concentration study. DLLME extracts of the cells treated with 0 (control), 50, 150, 300 and 1000 nM ergosterol, and incubated for 
18 h. The scores plots shows the clustering/separation of different treatments with little variation within each group. 
doi:1 0.1 371 /journal.pone.0087846.g001 
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Figure 2. Principal component analysis (PCA) loadings plots of 'H NMR spectral data. DLLME extracts were prepared from cells treated 
with 0 (control), 50, 150, 300 and 1000 nM ergosterol, and incubated for 18 h. The loadings plot shows the discriminating variables responsible for 
sample clustering depicted by the scores plots in Fig. 1 (B and C). (A) PCI vs. PC2 and (B) PC1 vs. PC3. 
doi:1 0.1 371 /journal.pone.0087846.g002 
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Figure 3. GC-FID chromatograms of DLLME extracts of ergosterol-treated (T) tobacco cells. (A) Chromatograms showing the differences 
between the control/non-treated cells (0 nM) and ergosterol-treated cells (50-1000 nM) incubated for 18 h. From bottom to top: 0 (control), 50, 150, 
300 and 1000 nM. (B) Chromatograms showing time-dependent variations of extracts of tobacco cells treated with 300 nM ergosterol and incubated 
for different time periods (0, 6, 12, 18 and 24 h T). The bottom chromatogram is a non-treated (NT) sample incubated for 24 h (24 h NT). 
doi:1 0.1 371 /journal.pone.0087846.g003 



treated (controls) and the 1000 nM ergosterol-treated samples are 
significantly separated from the other samples (i.e. different in 
metabolite composition). 

For the time study data, a three-component model was 
calculated and explained 89.6% variation [R 2 X(cum) of 0.896 
and Qf (cum) of 0.814]. The first two components (PCI and PC 2) 
explained 83.3% of the variance and were used to construct a 
scores plot [R 2 X{cum) of 0.833 and Qf(cum) of 0.742; with 95% 
confidence], showing samples differentially clustered into different 
groups (Figure 4B). Although treatment-related clusters formed, 
the 0 h treated samples were separated from the rest of the 
samples, while the non-treated samples (incubated for 24 h to 
compensate for aging of cells) were seen to group with the 6- and 
12 h-treated samples, and the 18 h- and 24 h-incubated samples 
were clustered together, indicating that the biggest changes 
occurred in the period 12-18 h, and that the response was 
essentially complete at 18 h. 

The GC-FID results thus reveal that ergosterol induces changes 
in DLLME-extractable secondary metabolites of tobacco cells. 
These changes include variation in the levels of the constitutively 



expressed metabolites and production of new metabolites, as 
demonstrated in Figures 3-A and -B. Differential changes in the 
intracellular metabolite profiles, explained by PCA (Figure 4), 
reflect the cells' response to the perception of ergosterol. 

Although the multi-linear programmed-temperature conditions 
in GC analyses improves the separation of the components in a 
sample [65,66], the ultracomplexity/multi-dimensionality of the 
extracted plant metabolites makes baseline separation difficult 
[3,67]. Overlapping peaks or co-eluting compounds that exist in 
the obtained chromatogram, can contribute to the observed 
variations. Furthermore, technical and instrumental variation can 
also contribute to minor variations in the chromatograms and thus 
influence the MVDA output [66]. Thus, to add a different 
dimension, DLLME samples were further analysed with a 
comprehensive GCxGC-TOF-MS technique. 

Two-Dimensional Gas Chromatography (GCxGC-TOF- 
MS) Analyses 

The advantage of the GCxGC-TOF-MS technology is its 
multi-dimensional separation, which is based on the orthogonality 
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Figure 4. PCA scores plots of the GC-FID concentration- and time study data. (A) DLLME extracts of cells treated with 0 (control), 50, 1 50, 
300 and 1000 nM ergosterol, and incubated for 18 h. The scores plot shows the clustering of different treatments (0-1000 nM) with little variation 
within each group. The 50-300 nM clusters appear to group together. The 0 nM (control) samples are clearly separated from the treated samples. (B) 
DLLME extracts of cells treated with 300 nM ergosterol and incubated for different time periods (0-24 h T) and non-treated samples incubated for 
24 h (24 h NT). The scores plot shows the clustering of the different treatments, with the 18 h- and 24 h-treated samples being significantly 
separated from the other samples. 
doi:1 0.1 371 /journal.pone.0087846.g004 



of two columns that are used in this technique [42]. GCxGC- 
TOF-MS permitted an improved chromatographic separation of 
the DLLME samples for both the time and concentration studies, 
compared to the one-dimensional GC-FID analyses by providing 
an expanded separation space and therefore minimising possible 
co-elution of analytes or 1-D peak overlap. Figure 5 is a 
representative of the TIC chromatograms obtained from the 
GCxGC-TOF-MS analyses. Generally, each spot/dot in the 2D/ 
3D-TIC chromatogram represents a detected ion peak. The visual 
inspection of these TIC chromatograms shows clear differences (in 
terms of the number/intensity of detected ion peaks) between 
control (non-treated) and ergosterol-treated samples. 

Although the DLLME method selectively extracts semi- and 
non-polar compounds, the resulting analytes represent only part of 
the total metabolome. DLLME extracts are still multi-dimensional 
samples due to the inherent ultracomplexity of the plant 
metabolome. Due to differential and multiple decoration of a 
common skeleton (by methylation, prenylation, hydroxylation, 
conjugation, or acylation reactions), plant secondary metabolites 
are structurally highly diverse, forming a highly complex spectrum 
of compounds of different size, solubility, volatility, polarity, 
quantity and stability [4,15,63,68]. 

Sample dimensionality is a measure of sample complexity and 
strongly influences component resolution in relationship to the 
dimensionality of the separation system [69] . Sample dimension- 
ality, denoted as s, is an intrinsic property of analytical samples 
(other than the number of constituents) that determines their 
amenability to separation systems. This parameter s is defined as 
the number of independent variables that must be specified to 
identify the components of the sample, and it arises mostiy from 
the physico-chemical differences of the sample constituents 
[69,70]. Semi- and non-polar plant metabolites present in the 
DLLME extracts would come from different chemical families, 
with diverse decoration, forming a multi-dimensional sample. 
Multi-dimensional and/or hyphenated separation systems would 
ideally provide a better separation of such complex samples 
compared to one-dimensional system, such as ID GC-FID [69- 
72]. 

The use of a GC xGC-TOF-MS in this study thus provided an 
enhanced peak capacity for gas chromatographic analyses of the 



DLLME extracts, and also revealed the internal structure of the 
sample constituents as represented by the 2D pattern of peak 
distribution in TIC chromatograms (Figure 5). The peak pattern 
generated from a GC x GC analysis can be used as defining the 2D 
signature of a sample's components [73] and can be useful for a 
metabolic fingerprinting study. The visual inspection of the 
obtained chromatograms was sufficient to provide a level of 
analytical certainty regarding the ergosterol-induced variations, 
which were already observed from ID GC-FID data. The MS- 
based annotation of compounds was obtained from the GC xGC- 
TOF-MS analyses of DLLME samples and is presented below. 

Generally, only relatively small molecules of low polarity can be 
directly analysed by GC; hence the ergosterol-induced variations 
depicted by the GC analyses is, to a certain degree, limited by the 
GC capability. Liquid chromatography, on the other hand, 
permits separation of compounds of a wide range of polarity 
and molecular weight [64,74,75]. Thus, to obtain an even more 
comprehensive view of the effect of ergosterol on the metabolome, 
LC coupled to high definition mass spectrometry (UHPLC- 
HDMS) technology was subsequently employed as analytical 
technique. Aspects of the UHPLC-HDMS analyses of DLLME 
samples was previously reported, and showed that ergosterol 
triggered differential changes in the non— /semi-polar metabo- 
lome of tobacco cells, with the de novo biosynthesis of five bicyclic 
sesquiterpenoids phytoalexins: capsidiol, lubimin, rishitin, solave- 
tivone and phytuberin [33]. The LC analyses of methanol 
extracted (ME) samples in this report also indicated differential 
metabolic changes related to ergosterol treatment (Figures S1-S3). 

Considering the analytical range of GC and LC, the results 
from both GC and LC analyses evidently suggest that the defence 
response induced by ergosterol in tobacco cells comprises a 
metabolic reprogramming that involves a wide range of intracel- 
lular secondary metabolites (from thermally stable non-polar 
metabolites to the thermally labile, semi-polar metabolites, with 
different molecular weights and chemical structures). Due to the 
ultracomplexity of the plant metabolome and the inherent 
limitations of the existing technologies, no single analytical 
approach is entirely competent in fully covering the whole 
metabolome or a specific metabolic range within a given biological 
sample [3-5,15]. Thus, the combination of ID ] H NMR, GC- 
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Figure 5. Representative 2D and 3D TIC chromatograms from the GC GC-TOF-MS analyses. The representative 2D/3D TIC 
chromatograms of the DLLME samples show differences in terms of the number of ion peaks detected and the variation in peak intensities between 
non-treated and 300 nM ergosterol-treated samples, incubated for a period of 18 h. 
doi:1 0.1 371 /journal.pone.0087846.g005 



FID, GC xGC-TOF-MS and UHPLC-MS technologies and two 
different extraction methods (a selective DLLME and general ME) 
provided a relatively comprehensive snapshot of the effect of 
ergosterol on the metabolome of tobacco cells. In order to 
elucidate clear, specific biological insights into the effect of 
ergosterol on plant metabolism, the study further translated the 
observed ergosterol-induced metabolic changes (demonstrated by 
visual inspection of chromatograms/NMR-spectra and PCA- 
explained variations), into the annotation of known metabolites. 

Metabolite Annotation 

Metabolite annotation or tentative/putative identification is one 
of the essential steps in metabolomic analyses, since biochemical 
interpretation of metabolomic data relies on the chemically 
annotated molecules from spectrometric/spectroscopic signals 
[12,14]. However, one of the difficulties that arise in metabolomic 
studies is the identification of the de nooo-induced compounds due 
to both the very restricted amounts and the high complexity of the 
biological extracts [76]. Furthermore, in the case of plant : 
pathogen interaction studies, the induced metabolites have varying 
rates of accumulation, and most of these metabolites may not be 
stable and can undergo degradation either by the plant or in the 
extraction process. In this regard, plants contain enzymes that 
degrade or convert antimicrobial compounds, returning their 
levels to pre-infection or pre-stress concentrations after the 
infection or stress has been contained or accommodated [77]. 



Thus, mass spectral data (from GC-EI-MS, GCx GC-TOF-MS 
and UHPLC-HDMS) were used for identification of biomarkers 
and de novo metabolites related to the observed ergosterol-induced 
changes in the tobacco cells. 

The mass spectral data provide a pattern that is often 
compound specific. However, the degree of certainty in elucidat- 
ing structural and chemical identity of an MS-detected analyte 
relies on the efficiency and accuracy of the three principal 
processes of the MS (ionisation, ml z analysis and ion detection) 
and on appropriate algorithms. The usage of parallel analytical 
platforms can provide additional information or confirmation for a 
putatively identified metabolite [1 1,12,78,79]. Here, accurate mass 
measurement, high mass spectrum match factors, accuracy in 
calculating elemental composition, usage of some MS E informa- 
tion (where available) and confirmation from parallel analytical 
platforms provided the certitude in metabolite assignments/ 
annotations, thus indicated as 'tentative' or metabolite identifica- 
tion (Ml)-level 2 [80]. 

GC-EI-MS and GCxGC-TOF-MS Based Annotation of Non- 
polar, Volatile Compounds 

For GC-EI-MS analysis, the search results against the NIST 
and Wiley mass spectra libraries were evaluated based on the 
probability factor, a match factor provided by the search 
algorithms that describes how unique a spectrum is, compared 
to all other spectra in the library [81-83]. Only spectra with 
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Table 1. Some of the annotated/tentatively identified (Ml-level 2) metabolites (with VIP scores >1.0) from ergosterol-treated 
tobacco cells. 



# 


m/z 


Rt (s) 


MF 


MW 
(ave) 


Compound ID name 


Platform* 


Metabolite 
category 


1 


265.1545 


599.4 


C,s H 24 0 


220.35 


Deoxy-capsidiol 


LC-MS 


Terpenoid 


2 


236 


990; 1.79 


C 15 0 24 0 2 


236 


Capsidiol 


LC-MS; 2D GC-MS 


Terpenoid** 


3 


204 


985; 1.65 


C15H24 


204 


epi-Aristolochene 


LC-MS; 2D GC-MS 


Terpenoid 


4 


222 


1009; 1.81 


C 15 H 26 0 


222 


Farnesol 


1 D GC-MS;2D GC-MS 


Terpenoid 


5 


155.1427 


745.8 


*-10 n 18 1 -' 


1 54.249 


Goran iol 


LC-MS 


Terpenoid 


5 


449.1 883 


616.2 


l -20 rl 36 l -'7 r 2 


450.443 


Gerany Igeranyl diphosphate 


LC-MS 


Terpenoid 


7 


291.1369 


595.8 




254.365 


H yd roxy-d i hyd ro I u b i m i n 


LC-MS 


Terpenoid 


8 


297 1 44 


1 85 4 


r r n nJ n, 

, -15 n 24 l - , 3 


252 349 


M\//~irr\Y\/li 1 him i n 

n y U 1 L/ Ay 1 LI )J \ II III 1 


LC-MS 


1 tri l-'tr 1 \\j 1 u 


9 


311 1 1 66 


642 6 


y ~^5 ^t 24'~'2 


236 349 


1 1 1 hfmin 

LU Ull 1 111 1 


LC-MS 


Tpi rnp n n i H 
1 d [Jt: 1 1 kj 1 Kx 


1 0 


261 1245 


1 75 2 


, -14 n 22 l -'2 


222 323 


Richitin 

Rl jl 1 1 LI 1 1 


I C-MS- 1D GC-MS 
i_v_ ivi i vj , zu uv_ ivij) 


1 L. 1 1 1 \J 1 u 


1 1 


309.144 


444 


*~^5 ^l 22 l ^ J 


218.334 


Solavetivone 


LC-MS; 2D GC-MS 


Terpenoid 


1 2 


389.1 81 


738 




390.425 


( — ) - 1 1 -Hydroxy-9, 1 0 - dihydrojasmonic acid 11 -beta-D- 
glucoside 


LC-MS 


Signaling 


13 


337.09 


570.6 


Q5H20O5 


280.316 


8'-Hydroxyabscisate 


LC-MS 


Signaling 


14 


295.1142 


774.6 


C,2H, 8 0 4 


226.268 


epi-4'-Hydroxyjasmonic acid 


LC-MS 


Signaling 


15 


211.1321 


597.6 


C,2H 18 0 3 


210.269 


Jasmonic acid 


LC-MS 


Signaling 


16 


299.078 


162 


C 13 H, 6 0 8 


300.261 


Salicylate 2-O-beta-D-glucoside 


LC-MS 


Signaling 


17 


385.116 


401.4 


C17H22O10 


386.35 


1 -O-Sinapoyl-beta-D-glucose 


LC-MS 


Phenylprop 


18 


195.066 


604.2 


C10H12O4 


196.199 


2'.6'-Dimethoxy-4'-hydroxyacetophenone 


LC-MS 


Phenylprop 


19 


319.083 


649.2 


C,6H, 6 0 7 


320.294 


4-Coumaroylshikimate 


LC-MS 


Phenylprop 


20 


177.055 


601.8 


C 10 H 8 O3 


176.168 


4-Methylumbelliferone 


LC-MS 


Phenylprop 


21 


351.071 


444.6 


C 16 H, 6 0 9 


352.292 


4-Methylumbelliferone glucuronide 


LC-MS 


Phenylprop 


22 


337.0925 


570.6 




338.309 


4-p-Coumaroylquinic acid 


LC-MS 


Phenylprop 


23 


341.088 


349.2 


CisHisOg 


342.298 


Caffeic acid 3-glucoside; 1-Caffeoyl-beta-D-glucose 


LC-MS 


Phenylprop 


24 


355.102 


459.6 




354.308 


Caffeoylquinate 


LC-MS 


Phenylprop 


25 


335.077 


559.8 


C 16 H l6 0 8 


336.293 


Caffeoylshikimate 


LC-MS 


Phenylprop 


26 


312.123 


780 


C 17 H l7 N0 2 


267.322 


Cinnamoyltyramine 


LC-MS 


Phenylprop 


27 


181.086 


492 


C 10 H l2 O 3 


180.2 


Coniferyl alcohol 


LC-MS 


Phenylprop 


28 


147.044 


570 


C 9 H 6 0 2 


146.142 


Coumarin 


LC-MS 


Phenylprop 


29 


515.117 


748.8 


^25H 2 40l2 


516.45 


Dicaffeoylquinic acid 


LC-MS 


Phenylprop 


30 


469.1107 


748.8 


C 21 H 22 0 8 


402.394 


Flavanone 7-O-beta-D-glucoside 


LC-MS 


Phenylprop 


31 


193.051 


699.6 


C10H10O4 


194.184 


Isoferulic acid 


LC-MS 


Phenylprop 


32 


251.151 


445.2 


C 12 H 18 N 4 0 2 


250.297 


N-caffeoylputrescine 


LC-MS 


Phenylprop 


33 


391.1 


376.8 


C17H20O9 


368.335 


O-Feruloylquinate 


LC-MS 


Phenylprop 


34 


191.057 


570.6 


C 7 H 12 0 6 


192.166 


Quinate 


LC-MS 


Phenylprop 


35 


193.049 


384 


CioH 8 0 4 


192.168 


Scopoletin 


LC-MS 


Phenylprop 


36 


173.046 


444.6 


C 7 H 10 O 5 


174.151 


Shikimate 


LC-MS 


Phenylprop 


37 


415.124 


522.6 


C17H22O9 


370.351 


Sinapaldehyde glucoside 


LC-MS 


Phenylprop 


38 


225.077 


492 


C,,H,20 5 


224.209 


Sinapate 


LC-MS 


Phenylprop 


39 


342.134 


796.8 


C 19 H 21 N0 5 


343.373 


Sinapoyltyramine 


LC-MS 


Phenylprop 


40 


373.148 


784.2 


C^H^Og 


372.367 


Syringin 


LC-MS 


Phenylprop 


41 


355.103 


642.6 


C15H18O7 


310.299 


trans-Cinnamoyl beta-D-glucoside 


LC-MS 


Phenylprop 


42 


195.065 


586.2 


C 8 H 14 0 4 


174.194 


Octadecanoic acid/suberic acid 


LC-MS; 2D GC-MS 


Fatty acid 



•Identification methods/platforms: ID GC-MS = GC-EI-MS; 2D GC-MS = GC xGC-TOF-MS; LC-MS = U H PLC-HDMS. 

**The 1 D proton NMR analyses of DLLME samples indicated differential changes in spectral region 3-5.8 ppm, which could be resonance from sesquiterpenoids (rishitin 
and capsidiol). 

doi:1 0.1 371 /journal.pone.0087846.t001 



probability hits in the range from 650 to 900 or above, were annotated from the GC-EI-MS analyses of DLLME samples 
considered. Six metabolites (terpenoids and fatty acids) were thus (Table 1 and Table SI). 



PLOS ONE | www.plosone.org 



10 



January 2014 | Volume 9 | Issue 1 | e87846 



Metabolomics of Ergosterol-lnduced Plant Responses 



For GCxGC-TOF-MS analyses, the Chroma-TOF software 
was used for automated data analysis from preprocessing to the 
compound identification step. The software, with its peak finding 
and deconvolution algorithms features, detects and extracts peaks 
with their "pure" mass spectra. ChromaTOF software identified 
compounds by comparing the extracted mass spectra to the NIST 
mass spectral library. A peak table was then generated from the 
library search and includes compound names, empirical formulae, 
Rt (in both dimensions), peak area, CAS numbers, mass spectral 
match factors and formula weights. The search results (the 
generated peak table) were evaluated based on a similarity match 
factor. The latter describes how well the library hit matches the 
peak when using all masses, and the values above 800 indicate the 
best fit [83]. In this study, the annotated compounds (similarity 
factor >800) from GC xGC-TOF-MS analyses include terpenoids 
and phytosterols (Table 1 and Table SI). 

UHPLC-MS Based Annotation of Mixed Polarity 
Compounds: In-depth Analysis of Metabolomic Data 

The variation evidenced by visual inspection of UHPLC-MS 
BPI chromatograms and explained by PC A modelling of UHPLC- 
MS data revealed differential changes in intra-cellular metabolite 
profiles, reflecting thus the response of tobacco cells to ergosterol 
treatment as a reprogramming of the metabolome [33]. In order 
to complement the view provided by PCA modelling and closely 
detect discrimination between ergosterol-treated samples and non- 
treated controls, a supervised model, OPLS-DA was used (for both 
ESI positive and negative data, from both DLLME and methanol, 
ME, samples). OPLS-DA differs from PCA by addition of 
grouping variables that indicate in which class the samples belong. 
Where PCA modelling is a descriptive method, OPLS-DA method 
is an explicative/predictive analysis. The latter facilitates the 
identification of the metabolite ions responsible for the discrim- 
ination between groups. OPLS-DA models data according to a 
priori class information (such as treated vs. non-treated) assigned to 
samples before the analysis; and it is a suitable tool to extract 
information on changes/ differences in the molecular composition 
of samples under study [14]. OPLS-DA modelling of control and 
300 nM ergosterol-treated (incubated for 18 h) samples was 
performed, to separate multivariate relationships into predictive 
variation (related to ergosterol treatment) and orthogonal variation 
(unrelated to ergosterol treatment). Analysis of variance testing of 
cross-validated predictive residuals (CV-ANOVA), a diagnostic 
tool, was used to assess the reliability of the obtained OPLS 
models. A /j-value that is lower than 0.05 normally indicates a 
significant model [84,85] . 

For DLLME samples, both ESI positive and negative data, the 
calculated OPLS-DA models were highly significant with a Rvalue 
of 0.0006 (ESI positive) and 0.005 (ESI negative), and comprised 
one predictive and three orthogonal components, with a 
reasonable fit to the data (R 2 X >0.714, R 2 Y = 0.99 and 
O_ 2 >0.817). For ESI positive data, 27% of variation in the 
computed OPLS-DA model was due to ergosterol treatment 
whereas for ESI negative data, only 22% of variation was related 
to ergosterol treatment. In both cases, about 50% of variation in 
the model was orthogonal variation. For ME samples, ESI positive 
data, the calculated OPLS-DA model was highly significant with a 
/(-value of 8.8x10 , and comprised one predictive and one 
orthogonal component, with a reasonable fit to the data 
(R 2 X=0.62, R 2 Y = 0.996 and Q. 2 = 0.985). 55% of variation in 
the computed OPLS-DA model was due to ergosterol treatment 
and 7% was orthogonal variation. For ME samples, ESI negative 
data, the calculated OPLS-DA model was highly significant with a 
Rvalue of 9.2x10 1 , and comprised one predictive and two 



orthogonal component, with a reasonable fit to the data 
(R 2 X= 0.645, R 2 Y= 0.993 and Q 2 = 0.955). 23% of variation in 
the computed OPLS-DA model was due to ergosterol treatment 
and 4 1 % was orthogonal variation. 

The evaluation of OPLS-DA loading S-plots (Figure 6A) 
permitted the extraction of statistically and potentially biochem- 
ically significant mass ions (metabolites/bio-markers) in the 
samples. These mass ions, related to ergosterol treatment, were 
selected based on their contribution to the model (x-axis, modelled 
covariation) and reliability (y-axis, modelled correlation). An ideal 
statistically significant variable (mass ion, in this case) has high 
covariation (magnitude) and high reliability, i.e. smaller risk for 
spurious correlations [48,54]: \p [1] | >0.05 and |/>(corr)| ^0.5, 
respectively, in this study. Furthermore, the significance of the 
selected mass ions (variables) in OPLS-DA models was also 
assessed using the variable importance in projection (VIP) 
(Figure 6B). The VIP plot summarises the importance of the 
variables both to explain X and to correlate to Y. The VIP score is 
a critically important check on the selection of significant ions/ 
variables in a complex data set from metabolomics. The higher the 
VIP value (exceeding 1.0) the more significant is the ion/variable 
in the complex analysis in comparing difference between two or 
more groups [86-88]. These mass ions are presented in Table 1 
(VIP>1.0) and Table SI (VIP<1.0), together with their annota- 
tions/putative identities. 

The selected mass ions from S- and VIP plots were annotated 
using the Taverna workbench (www.taverna.org.uk) for PUTME- 
DID_LCMS Metabolite ID Workflows. The Taverna workflows 
allow for integrated, automated and high-throughput annotation 
and tentative metabolite identification from ESI-LC-MS meta- 
bolic data. The workflows consist of correlation analysis, metabolic 
feature annotation and metabolite annotation [89]. 

Ergosterol Induces Metabolic Reprogramming That 
Involves Terpenoid and Phenylpropanoid Pathways in 
Tobacco, Leading to the Biosynthesis of Defence-Related 
Secondary Metabolites 

For metabolomics to be meaningful and successful, raw 
analytical data should be converted to structurally elucidated 
compounds/ metabolites in order to provide biological knowledge 
about the system under investigation [12]. The visual inspection of 
chromatograms/NMR-spectra and PCA-modelling evidently sug- 
gest that the fungal derived MAMP, ergosterol, induced a 
metabolic reprogramming in tobacco cells that involves a wide 
range of intracellular secondary metabolites (from thermally stable 
non-polar metabolites to the thermally labile, semi-polar metab- 
olites, with different molecular weights and chemical structures). 
Furthermore, translating the observed ergosterol-induced meta- 
bolic changes into structurally elucidated metabolites demonstrat- 
ed that this cellular reprogramming involves a number of 
metabolic pathways, such as the terpenoid- and phenylpropanoid 
pathways and their branches. Although no fluxomics analyses 
were conducted in this study, the annotation of secondary 
metabolites that are positively or negatively correlated with the 
treatment, sufficiendy suggests dynamic intra- and interrelations 
among different metabolic pathways (Table 1; Table SI; Figure 7). 

The ergosterol-induced components of the terpenoid pathway 
include the bicyclic sesquiterpenoids (capsidiol, lubimin, phytu- 
berin, rishitin and solavetivone) and their biosynthetic precursor, 
aristolochene; sesquiterpenoid derivatives (such as deoxy-capsidiol 
and hydroxylubimin); farnesol, an acyclic sesquiterpene; other 
precursors in the terpenoids pathway (geranyl diphosphate, GPP) 
as well as phytosterols (campesterol, sitosterol and stigmasterol) 
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Figure 6. OPLS-DA modelling of data from ergosterol-treated cells. (A) S-plot of an OPLS-DA model of control vs. 300 nM ergosterol-treated 
sample extracts (ME/UHPLC-MS data, ESI positive). The x-axis is the modelled covariation (variable magnitude) and the y-axis is the loading vector of 
the predictive component (modelled correlation/reliability). The mass ions in the upper right quadrant of S-plot are positively related to the 
ergosterol treatment (these include ions m/z 265.15, 177.051, 147.044 and 181.086), while those in the lower left quadrant are negatively related to 
the treatment. (B) A variable importance in projection (VIP) plot for an OPLS-DA model of 300 nM-treated samples. The VIP plot indicates for instance 
that the mass ions m/z 256.15, 177.051, 147.044, 295.164 and 181.086 (identified as deoxy-capsidol, methyllumbeliferone, coumarin, epi- 
hydroxyjasmonic acid and coniferyl alcohol, respectively, Table 1) were accountable for the significant separation in the model as their VIP scores 
were significantly greater than 1.0. 
doi:1 0.1 371 /journal.pone.0087846.g006 



and their precursor, squalene (Table 1; Table SI; Figure 7). These 
bicyclic sesquiterpenoid metabolites are found in plants within the 
Solanaceae, and are correlated with the defence responses to 
invading pathogens [90,91]. Furthermore, farnesol could be 



conceived of as being involved in regulating the rate of formation 
of other defence-related terpenoids [92]. 

The terpenoids are generally synthesized from the isomeric 5- 
carbon building block molecules, isopentenyl diphosphate (IPP) 



PLOS ONE | www.plosone.org 



12 



January 2014 | Volume 9 | Issue 1 | e87846 



Metabolomics of Ergosterol-lnduced Plant Responses 



^RIMARY METABOUS^ 




Monoterpenes 



GPP 



Sterols: 
Campesterol 
Sitosterol 
Stigmasterol 



2 x IPP 
. Squalene < FPP 



DMAPP 

3 x IPP 

GGPP 



1° 



1! 



)iterpertes 
Aristolochene ,:, 



Debneyol 



: '•..Jetraterpenes 

• Capsidiol 



Solavetivone ' 

^ v 

▼ Phytuberin 
Lubimin ^^•Hydroxylubimin Phytuberol 

Rishitin 



Dehydroquinate " 



Dehydroshikimate 
Shikimate 



Chorismate 
Phenylalanine 



Cinnamic acid (CA) 
H/M CAs 

H/M CA-CoA esters 



quinate 



COUMARINS 
4 



B-D-glu-CA 



ACETOPHENONES 

0 



p-coumarate 
Caffeate 
Ferulate 
Sinapate 



H/M cinnamoyl alcohols 
H/M cinnamoyl aldehydes 



p-coumaryl alcohol 
p-coumaraldehyde 
coniferyl alcohol 
coniferaldehyde 
sinapyl alcohol 
sinapaldehyde 



V 



BENZOIC ACIDS 



.r>.FLAVONOIDS 



'4 



CONJUGATES 

agmatine 

tyramine 

p-coumaroyl quinate 
p-coumaroyl shikimate 
caffeoyl quinate 
caffeoyl shikimate 



..> LIGNIN 



Figure 7. The ergosterol-induced metabolic reprogramming in N. tabacum cells. The fungal derived MAMP, ergosterol, induced an intra- 
cellular metabolic reprogramming that involves terpenoid and phenylpropanoid pathways and their branches, indicated on this diagram. IPP, 
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methoxylated cinnamic acids. Metabolites shown in blue have been tentatively identified in this study (Table 1 and Table SI). 
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and dimethylallyl diphosphate (DMAPP), the products of two 
independent pathways in plants; the mevalonate (MVA) pathway 
operating in the cytosol and the 2-C-methyl erythritol 4-phosphate 
(MEP) pathway in plastids [93-96] . A series of enzyme-catalysed 
condensation reactions of IPP and DMAPP molecules leads to the 
biosynthesis of farnesyl diphosphate (FPP), a 15-carbon molecule. 
FPP is further converted into a sesquiterpene, 5-«/«'-aristolochene, 
by the action of a sesquiterpene synthase. Sesquiterpenoids are 
subsequently enzymatically derived from this common precursor 
[94,97-99]. 

The five sesquiterpenoids annotated in this study (also in [33]) 
have been reported to accumulate in plant cell suspension cultures 
or tissues challenged by fungal elicitors or by viral infection, 
thereby providing an antimicrobial and fungitoxic environment 
[90,100,101]. These biotic stress-induced sesquiterpenoids are 
known phytoalexins [91,98,101], referring to host-synthesised low 
molecular weight compounds with protective properties of which 
de novo biosynthesis and accumulation are induced in plants 
following biotic or abiotic stress. The phytoalexin response is an 
effective part of the multi-component stress resistance mechanism 
in plants [101-103]. 

Furthermore, a relative quantitative analysis of selected 
sesquiterpenoids (rishitin, solavetivone and phytuberin) (Figure 8- 
A, -B and -C, respectively) indicated that accumulation of these 
phytoalexins was significantly detectable from 12 h post-treat- 



ment, reached a maximum at 18 h (/)<0.0001) and at 24 h their 
levels started to decrease. It is possible that from 24 h, these 
induced metabolites are either secreted to the extracellular milieu 
or are degraded. After the infection or stress has been contained or 
accommodated, the plant normally degrades antimicrobial com- 
pounds in enzyme-catalysed reactions, returning their levels to 
pre-infection or pre-stress concentrations [77]. 

Some signaling molecules and -derivatives/ conjugates were also 
tentatively identified. Plant hormone signaling pathways related to 
defence are not isolated but rather interconnected through 
complex regulatory networks. The type of interactions and plant 
responses to stresses vary depending on the pathosystem as well as 
the time, quantity and the tissue where the hormones are 
produced. To understand how plants coordinate multiple 
hormonal components in response to various developmental and 
environmental cues is a major challenge [104]. The annotated 
hormones include jasmonic acid (JA) and its hydroxylated 
derivative; salicylate (SA - a known defence signal) glucoside; 
methylsalicylate (MeSA); as well as abscisate (ABA) and hydro- 
xyabscisate, an ABA degradation product (Table 1, Table SI). 
Studies have shown that the activity of ABA-8' -hydroxylase is up- 
regulated by ABA accumulation and stress [105]. ABA is an 
isoprenoid plant hormone (C 1 5 -molecule) synthesised in the 
plastidial MEP pathway [106,107], unlike other structurally 
related C 15 -terpenoids (sesquiterpenoids) that are synthesised from 
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Figure 8. The quantitative, time-dependent, responses of some of the annotated compounds from different metabolite classes: (A- 

C) sesquiterpenoid phytoalexins (rishitin, solavetivone and phytuberin), (D-E) sterols (squalene and stigmasterol) and (F-H) phenolics (shikimate, 
caffeoyl shikimate and coniferyl alcohol). This relative quantification is based on their respective peak areas estimated from peak integration. The 
quantitative estimation shows a time-dependent biosynthesis of sesquiterpenoids, sterols and phenylpropanoids from cells treated with 300 nM 
ergosterol, reaching a maximum at 18 h post-treatment (N-ANOVA t-test, p<0.0001). 
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the cytosolic MVA pathway. ABA is vitally involved in the 
regulation of various aspects of plant growth and development. 
Furthermore, ABA regulates plant defence responses negatively or 
positively, depending on the type of plant : pathogen interaction or 
the nature of abiotic stress. The role of ABA in plant resistance 
appears to be complex and the exact molecular mechanism of this 
phytohormone in defence responses is still unclear [104,108-1 13]. 
The findings from a study conducted on Nicotiana plumbagintfolia 
suggested that even though ABA may not be required per se for the 
biosynthesis of capsidiol (a sesquiterpenoid phytoalexin), this 
phytohormone is vitally implicated in a stress-response checkpoint 
to fine-tune the amplification of capsidiol synthesis in challenged 
plants [1 14]. Azelaic acid, a mobile molecule that can confer local 
and systemic resistance via priming of plants to accumulate SA 
upon infection [115], was also found in extracts from ergosterol- 
treated cells (Table SI). 

The MS-based annotation of treatment-correlated metabolites 
indicates further that the ergosterol-induced metabolic reprogram- 
ming involved also the phenylpropanoid pathway and its branches 
(Figure 7; Table 1; Table SI). Coumarins, acetophenones and 
benzoic acids as well as a range of lignin precursor molecules were 
found. These metabolites are known in the context of plant stress - 
and defence responses [68]. Of special interest is the different 
conjugates of cinnamic acid derivatives that were annotated as 
being positively correlated with the ergosterol-induced responses. 
A close examination of some of the annotated phenylpropanoid 
compounds indicated also a significant time-dependent accumu- 
lation of these compounds (p<0.0001) (Figures 8-F to -H). The 
hydroxylated and methoxylated cinnamoyl-tyramine conjugates 
have been reported to accumulate in response to priming with 
bacterial lipopolysaccharides (LPS) and pathogen inoculation. 
These compounds have antimicrobial activity and can also be 
cross-linked into the plant wall to strengthen it and can thus play a 
role in the resistance response of plants [20,116]. The results 
indicate that ergosterol affects a wide activation of the shikimate - 
and phenylpropanoid - linked defence pathways, the products of 
which may subsequently lead to an antimicrobial environment 
in vivo. 

The fatty acids (FAs, saturated and unsaturated) annotated from 
GC-MS and GC xGC-TOF-MS analyses (Table 1; Table SI), 
include n-hexadecanoic acid (16:0, palmitic acid), octadecanoic 
acid (18:0, stearic acid) and 9, 1 2-octadecadienoic acid (Z,Z) (18:2, 
linoleic acid). Recent studies have demonstrated that FAs and their 
breakdown products participate actively and directly in various 
modes of plant defences, such as remodelling cell membrane 
fluidity and modulating defence gene expression. Both C 16 and 
C 18 FAs are involved in regulating basal, effector-triggered, and 
systemic immunity in plants. It has been shown that oleic acid 
(18:1) and linoleic acid (18:2) play significant roles in the activation 
of NADPH oxidase, resulting in the production of reactive oxygen 
species (ROS), both components of the plant defence reaction 
[117-119]. 

Phytosterols (campesterol, stigmasterol and P-sitosterol) and 
their precursor, squalene were among the compounds annotated 
using GCxGC-TOF-MS (Table 1; Table SI). Relative quantifi- 
cation indicated that squalene accumulated significantly 1 2 h post- 
treatment (p<0.0001) and the level decreased as that of 
stigmasterol reached a maxium at 18 h (/><0.0001) (Figures 8-D 
and -E). This observed variation of the detected phytosterols' level 
suggests a time-dependent response to ergosterol treatment. 
Although the sterols are constitutive major components of the 
cell membrane, recent discoveries (in jV. benthamiana and Arabidopsis 
thaliana systems) have demonstrated that plants alter sterol 
biosynthesis (particularly stigmasterol and R-sitosterol) upon 



pathogen attack so as to restrict nutrient transfer from the cytosol 
to the apoplastic space [120]. Furthermore, studies have pointed 
out changes in expression levels of A. thaliana C24 methyltrans- 
ferases Atlg20330 and At5gl3710 genes following the treatment 
with LPS [121]. These genes encode enzymes that control carbon 
flux into sterol biosynthesis, influencing subsequently membrane 
permeability. A pathogen-induced alteration in the stigmasterol/ 
R-sitosterol ratio can influence the physicochemical properties of 
membrane microdomains and thereby modulate plant defence 
signaling [122]. 

Both sesquiterpenoids and sterols are biosynthesised via the 
terpenoid pathway. Studies have shown that the position of FPP in 
the terpenoid biosynthetic pathway is an important and potential 
regulatory branch point of sesquiterpenoid biosynthesis. Under 
normal conditions, the FPP is channelled toward the biosynthesis 
of sterol and prenyl-lipid moieties [94,99]. Moreover, the 
enhanced expression of sesquiterpene cyclase, a key enzyme in 
the biosynthesis of sesquiterpenoids in response to treatment with 
tobacco mosaic virus or cell wall fragments from Phytophthora 
species, was accompanied by the suppression of squalene synthase 
activity in tobacco and potato cells/tissues [32,94,123,124]. Thus, 
the detection of both sesquiterpenoids and phytosterols in response 
to ergosterol suggests a complex regulation of the terpenoid 
pathway, permitting the production of both the sesquiterpenoid 
phytoalexins (fungitoxicity) and the increase of the levels of 
phytosterols (to restrict cell membrane permeability). However, 
further quantitative analyses and molecular-based studies are 
needed to corroborate these observations and unfold the molecular 
mechanisms involved in the complexity of regulatory network 
structures. 

Conclusion 

Ergosterol acts as a MAMP molecule in tobacco and tomato 
plants, triggering a defence response characterised by the 
elicitation of the oxidative burst, the production of ROS [31,32], 
alkalinisation of the external milieu and mobilisation of cytosolic 
calcium [30,31,125,126]. Furthermore, ergosterol induces expres- 
sion of some genes encoding pathogenesis-related proteins (PR1, 
PR3, and PR5) and enzymes participating in the defence response 
such as phenylalanine-ammonia lyase and sesquiterpene cyclase 
[32]. This ergosterol-induced defence response, a cellular repro- 
gramming, is reflected also by dynamic and differential changes to 
the metabolome as demonstrated by our results. 

The presence of the five bicyclic sesquiterpenoids (phytuberin, 
solavetivone, capsidiol, lubimin and rishitin) in ergosterol-elicited 
tobacco cells and other annotated metabolites (abscisic acid, fatty 
acids, phytosterols and metabolites of shikimate-phenylpropanoid 
pathways) indicates that the changes in the metabolome, 
demonstrated by chromatographic/spectral analyses and PCA- 
explained variations, are associated with a defensive function 
('defensome') in response to elicitation by ergosterol as a MAMP 
molecule. Ergosterol induced a complex and dynamic activation of 
the terpenoid, shikimate-phenylpropanoid and FA pathways 
leading to de novo biosynthesis of sesquiterpenoids phytoalexins 
and alteration in signaling molecules (JA, SA, MeSA, ABA, azelaic 
acid), FAs, phenylpropanoid-metabolies and phytosterols involved 
in various modes of defence responses. 

This ergosterol-induced metabolic reprogramming is dynamic 
and largely complex as it involves changes in different metabolic 
pathways which points to interconnected alterations in metabolic 
networks that are functionally correlated. The groups of identified 
compounds serve not only as a base in the search of novel defence 
compounds, but also as signatory bio-markers for the character- 
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isation of the plants' defensive state. This study identified mainly 
sesquiterpenoids, phytosterols, phenylpropanoid-metabolites (and 
conjugates/derivatives) and fatty acids, hence further metabolic 
network investigations are needed to unfold these intra- and 
interrelationships among different metabolic pathways involved in 
the ergosterol-induced response and to provide insights into the 
molecular mechanisms tangled in the complexity of regulatory 
network structures. 

Supporting Information 

Figure SI ESI positive BPI MS chromatograms of 
methanol extracts separated by UHPLC. (A) Concentra- 
tion study. The chromatograms show treatment-related variations. 
ME samples of the cells treated with different ergosterol 
concentrations (from bottom to top: 0, 50, 150, 300 and 
1000 nM) and incubated for 18 h. (B) Time study. The 
chromatograms show some time-dependent treatment-related 
variations. ME samples of cells treated with 300 nM ergosterol 
and incubated for different time periods (from bottom to top: 0, 6, 
12, 18 and 24 h T) and the top chromatogram is a non-treated 
sample incubated for 24 h (24 h NT). Visual inspection of these 
compared mass chromatograms (in both sets: concentration- and 
time study) indicate differential changes in detected mass ions, 
reflecting ergosterol-induced changes (concentration- and time- 
dependent) in the (ME-) metabolite profiles of the tobacco cells. 
(TIF) 

Figure S2 PCA scores plot of methanol extracts sepa- 
rated by UHPLC-MS in ESI positive mode. (A) Concen- 
tration study. ME samples of the cells treated with 0 nM (control), 
50 nM, 150 nM, 300 nM, and 1000 nM ergosterol, and incubat- 
ed for 18 h. A PCA-four-component model was calculated and 
explained 77.3% of the variance [R 2 X(cum) of 0.772 and Qf(cum) 
of 0.522]. Using the first two principal components (PCI and PC2, 
explaining 69.5% of the variance, with 95% confidence level) for a 
scores plot, the tobacco samples were found to be differentially 
clustered into five groups corresponding to different ergosterol 
treatments (0-1000 nM) with no significant intra-group variation. 
The clusters corresponding to 150 nM, 300 nM and 1000 nM are 
seen to group close to each other. The control (non-treated) and 
50 nM samples are clustered close to each other and significandy 
separated from the rest of the treatments. This scores plot 
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evidenced a clear dosage-dependent variation, suggesting also a 
similarity and less variability in metabolite content of 0 nM and 
50 nM samples. (B) Time study. ME samples of the cells treated 
with 300 nM ergosterol, and incubated for incubated for different 
time periods (0 h T - 24 h T) and a non-treated sample incubated 
for 24 h (24 h NT). A two-component model was calculated and 
explained 68.7% of the variance [i^X(cum) of 0.688 and Qf(cum) 
of 0.580]. A PCA scores plot (95% confidence) was constructed 
using PCI and PC2, showing samples differentially clustered into 
different groups: the 6 h- incubated samples formed a separate 
cluster, and the rest of the samples clustered separately but closer 
to each other. The scores plot shows a time-dependent variation. 
(TIF) 

Figure S3 OPLS-DA model of non-treated vs. 300 nM 
ergosterol-treated samples (ME/UHPLC-MS data, ESI 
negative). The calculated model had one predictive and one 
orthogonal components, and explained 81.7% of the variation 
[68.7% treatment-related variation (predictive) and 13.0% or- 
thogonal variation]. (A) The loading S-plot of the calculated 
model. The x-axis is the modelled covariation (variable magnitude) 
and the j-axis is the loading vector of the predictive component 
(modelled correlation/ reliability). The mass ions in the upper right 
quadrant of S-plot are positively related to the ergosterol 
treatment, while those in the lower left quadrant are negatively 
related to the treatment. (B) A variable importance in projection 
(VIP) plot for an OPLS-DA model of 300 nM-treated samples. 
The VIP plot indicates for instance that the mass ions ml 'z 
351.069, 335.077 and 337.09 (identified as methylumbeliferone 
glucuronide, caffeoylshikimate and hydroxyabscisate, respectively, 
Table 1 and Table SI) were accountable for the significant 
separation in the model as their VIP scores were significandy 
greater than 1.0. 
(TIF) 

Table SI Some of the tentatively identified (M-level 2) 
metabolites (with VIP< 1 .0) from ergosterol-treated tobacco cells. 
(DOCX) 

Author Contributions 

Conceived and designed the experiments: IAD FT. Performed the 
experiments: FT PAS. Analyzed the data: FT PAS IAD. Contributed 
reagents/materials/analysis tools: IAD. Wrote the paper: FT LAP IAD. 



11. Moco S, Vcrvoort J, Bino R, Dcvos R (2007) Metabolomics technologies and 
metabolite identification. TrAC Trends Anal Chem 26: 855—866. 

12. Brown M, Dunn WB, Dobson P, Patel Y, Winder CL, et al. (2009) Mass 
spectrometry tools and metabolite-specific databases for molecular identifica- 
tion in metabolomics. The Analyst 134: 1322-1332. 

13. Bhalla R, Narasimhan K, Swarup S (2005) Metabolomics and its role in 
understanding cellular responses in plants. Plant Cell Rep 24: 562-571. 

14. Tugizimana F, Piatcr LA, Dubery IA (2013) Plant metabolomics: A new 
frontier in phytochemieal analysis. S AfrJ Sci 109 (5/6): 1-11. doi:10.1590/ 
sajs.20 13/20 120005. 

15. Dunn WB, Ellis DI (2005) Metabolomics: current analytical platforms and 
methodologies. Trends Anal Chem 24: 285-294. 

16. Fukusaki E, Kobayashi A (2005) Plant metabolomics: potential for practical 
operation. J Biosci Bioeng 100: 347—354. 

17. Weckwcrth W, Loureiro ME, Wenzel K, Fiehn O (2004) Differential metabolic 
networks unravel the effects of silent plant phenotypes. Proc Natl Acad Sci 101: 
7809-7814. 

18. Lopcz-Gresa MP, Maltese F, Belles JM, Conejero V, Kim HK, et al. (2010) 
Metabolic response of tomato leaves upon different plant-pathogen interac- 
tions. Phytoehem Anal 2 1 : 89-94. 

19. Kim HK, Choi YH, Vcrpoortc R (2011) NMR-bascd plant metabolomics: 
where do we stand, where do we go? Trends Biotechnol 29: 267-275. 

20. Ahferis KA, Jabaji S (2012) FT-ICR/MS and GC-EI/MS metabolomics 
networking unravels global potato sprout's responses to Rhizoctonia solani 
infection. PloS One 7: 1—13. 



PLOS ONE | www.plosone.org 



16 



January 2014 | Volume 9 | Issue 1 | e87846 



Metabolomics of Ergosterol-lnduced Plant Responses 



21. Ingham JL (1972) Phytoalcxins and other natural products as factors in plant 
disease resistance. Bot Rev 38: 343—4-24. 

22. Fischer R, Hain R (1994) Plant disease resistance resulting from the expression 
of foreign phytoalcxins. Gurr Opin Biotechnol 5: 125-130. 

23. Wink M (1988) Plant breeding: importance of plant secondary metabolites for 
protection against pathogens and herbivores. Theor Appl Genet 75: 225-233. 

24. Jackson AO, Taylor CB (1996) Plant-microbe interactions: life and death at the 
interface. Plant Cell 8: 1651-1668. 

25. Croteau R, Kutchan TM, Lewis NG (2000) Natural Products (secondary 
metabolites). In: Buchanan B, Gruisscm W., Jones R, editors. Biochemistry & 
Molecular Biology of Plants. USA: American Society of Plant Physiologists. 
1250-1318. 

26. Pictcrsc CMJ, Leon-Reyes A, Van dcr Ent S, Van Wees SGM (2009) 
Networking by small-molecule hormones in plant immunity. Nat Ghem Biol 5: 
308-316. 

27. Pritchard L, Birch P (2011) A systems biology perspective on plant-microbe 
interactions: biochemical and structural targets of pathogen effectors. Plant Sci 
180: 584-603. 

28. Sanabria NM, Goring D, Niirnberger T, Dubery IA (2008) Sclf/non-self 
perception and recognition mechanisms in plants; a comparison of self- 
incompatibility and innate immunity. New Phytol 178: 503-514. 

29. Sanabria NM, Huang J, Dubery IA (2009) Self/non-self perception in plants in 
innate immunity and defence. Self/Nonsclf: Immune Recognition and 
Signaling 1: 1—15. 

30. Kasparovsky T, Milat M, Humbert C, Blein J, Havel L, et al. (2003) Elicitation 
of tobacco cells with ergostcrol activates a signal pathway including 
mobilization of internal calcium. Plant Physiol Biochcm 41: 495-501. 

31. Kasparovsky T, Blein J-P, Mikes V (2004) Ergostcrol elicits oxidative burst in 
tobacco cells via phospholipase A2 and protein kinase C signal pathway. Plant 
Physiol Biochem 42: 429-435. 

32. Lochman J, Mikes V (2006) Ergostcrol treatment leads to the expression of a 
specific set of defence-related genes in tobacco. Plant Mol Biol 62: 43-51. 

33. Tugizimana F, Stcenkamp PA, Piatcr LA, Dubery IA (2012) Ergostcrol- 
induccd sesquitcrpenoid synthesis in tobacco cells. Molecules 17: 1698-1715. 

34. Van Gulik WM, Ten Hoopen HJG, Heijnen JJ (2001) The application of 
continuous culture for plant cell suspensions. Enzyme Microb Tcchnol 28: 
796-805. 

35. Ghoi H-K, Yoon J-H (2007) Metabolomic profiling of Vitis vinifcra cell 
suspension culture elicited with silver nitrate by *H NMR spectrometry and 
principal components analysis. Process Biochcm 42: 271—274. 

36. Florcs-Sanchcz IJ, Pec J, Fei J, Ghoi YH, Dusck J, et al. (2009) Elicitation 
studies in cell suspension cultures of Cannabis sativa L.J Biotechnol 143: 157- 
168. 

37. Gcrber IB, Dubery IA (2004) Protein phosphorylation in Nicotiana tabacum 
cells in response to perception of lipopolysaccharides from Burkholdcria 
cepacia. Phytochcmistry 65: 2957—2966. 

38. Sanabria NM, Dubery IA (2006) Differential display profiling of the Nicotiana 
response to LPS reveals elements of plant basal resistance. Biochem Biophys 
Res Commun 344: 1001-1007. 

39. Towill LE, Mazur P (1975) Studies on the reduction of 2,3,5-triphcnyltc- 
trazolium chloride as a viability assay for plant tissue cultures. Can J Bot 53: 
1097-1102. 

40. Weckwerth W, Fiehn O (2002) Can we discover novel pathways using 
metabolomic analysis? Curr Opin Biotechnol 13: 156-160. 

41. Xiao-huan Z, Qiu-hua WU, Mei-yue Z, Guo-hong XI, Zhi W (2009) 
Developments of dispersive liquid-liquid microcxtraction technique. 
Chinese J Anal Chcm 37: 161-168. 

42. Marriott P, Shellie R (2002) Principles and applications of comprehensive two- 
dimensional gas chromatography. Trends Anal Ghem 21: 573—583. 

43. Shellie RA, Wclthagcn W, ZrosdikovaJ, Spranger J, Ristow M, et al. (2005) 
Statistical methods for comparing comprehensive two-dimensional gas 
chromatography-timc -of- flight mass spectrometry results: metabolomic anal- 
ysis of mouse tissue extracts. J Chromatogr A 1086: 83-90. 

44. Kock MM, Van dcr Kloet FM, Kleemann R, Kooistra T, Verheij ER, et al. 
(2011) Semi-automated non-target processing in GCxGC-MS metabolomics 
analysis: applicability for biomedical studies. Metabolomics 7: 1-14. 

45. Canellas E, Nerin C, Moore R, Silcock P (2010) New UPLC coupled to mass 
spectrometry approaches for screening of non-volatile compounds as potential 
migrants from adhesives used in food packaging materials. Anal Chim Acta 
666: 62-69. 

46. Peters RJB, Bolck YJC, Rutgers P, Stolkcr AAM, Nielen MWF (2009) Multi- 
residue screening of veterinary drugs in egg, fish and meat using high- resolution 
liquid chromatography accurate mass time -of- flight mass spectrometry. 
J Chromatogr A 1216: 8206-8216. 

47. Allwood JW, Goodacre R (2010) An introduction to liquid chromatography- 
mass spectrometry instrumentation applied in plant metabolomic analyses. 
Phytochem Anal 21: 33-47. 

48. Trygg J, Holmes E, Lundstcdt T (2007) Chcmometrics in metabonomics. 
J Proteome Res 6: 469-479. 

49. Bylesjo M, Eriksson D, Sjodin A, Jansson S, Moritz T, ct al. (2007) Orthogonal 
projections to latent structures as a strategy for mieroarray data normalization. 
BMC Bioinfbrmatics 8: 1-10. 



50. Ni Y, Su M, Lin J, Wang X, Qiu Y, ct al. (2008) Metabolic profiling reveals 
disorder of amino acid metabolism in four brain regions from a rat model of 
chronic unpredictable mild stress. FEBS Letters 582: 2627-2636. 

51. Tulipani S, Llorach R, Jaurcgui O, Lopcz-Uriarte P, Garcia-Aloy M, et al. 
(2011) Metabolomics unveils urinary changes in subjects with metabolic 
syndrome following 12-week nut consumption. J Proteome Res 10: 5047-5058. 

52. Camacho D, Fucntc A De, Mendcs P (2005) The origin of correlations in 
metabolomics data. Metabolomics 1: 53—63. 

53. Van den Berg RA, Hoefsloot HCJ, Westerhuis JA, Smilde AK, Werf MJ Van 
Dcr (2006) Centering, scaling, and transformations: improving the biological 
information content of metabolomics data. BMC Genomics 7: 1-15. 

54. Wiklund S, Johansson E, Sjostrom L, Mellcrowicz EJ, Edlund U, et al. (2008) 
Visualization of GC/TOF-MS-based metabolomics data for identification of 
biochemically interesting compounds using OPLS class models. Anal Ghem 80: 
115-122. 

55. Choi YH, Kim HK, Hazekamp A, Erkclens C, Lefeber AWM, et al. (2004) 
Metabolomic differentiation of Cannabis sativa cultivars using ] H NMR 
spectroscopy and principal component analysis. J Nat Prod 67: 953-957. 

56. Grebenstein C, Choi YH, RongJ, De Jong TJ, Tamis WLM (201 1) Metabolic 
fingerprinting reveals differences between shoots of wild and cultivated carrot 
(Daucus carota L.) and suggests maternal inheritance or wild trait dominance in 
hybrids. Phytochemistry 72: 1341-1347. 

57. Kim HK, Choi YH, Verpoorte R (2010) NMR-based metabolomic analysis of 
plants. Nat Protoc 5: 536-549. 

58. Choi YH, Kim HK, linthorst HJM, Hollander JG, Lefeber AWM, et al. (2006) 
NMR metabolomics to revisit the tobacco mosaic virus infection in Nicotiana 
tabacum leaves. J Nat Prod 69: 742-748. 

59. Vichi M, Saporta G (2009) Clustering and disjoint principal component 
analysis. Gomput Stat Data An 53: 3194-3208. 

60. Zhao J, Davis LC, Verpoorte R (2005) Elicitor signal transduction leading to 
production of plant secondary metabolites. Biotechnol Adv 23: 283—333. 

61. Sanchcz-Sampcdro A, Kim HK, Ghoi YH, Verpoorte R, Corchete P (2007) 
Metabolomic alterations in elicitor treated Silybum marianum suspension 
cultures monitored by nuclear magnetic resonance spectroscopy. J Biotechnol 
130: 133-142. 

62. Vasconsuelo A, Boland R (2007) Molecular aspects of the early stages of 
elicitation of secondary metabolites in plants. Plant Sci 172: 861-875. 

63. Yazaki K, Sasaki K, Tsurumaru Y (2009) Prenylation of aromatic compounds, 
a key diversification of plant secondary metabolites. Phytochcmistry 70: 1739- 
1745. 

64. Kaal E, Janssen H-G (2008) Extending the molecular application range of gas 
chromatography. J Chromatogr A 1184: 43-60. 

65. Blumbcrg LM, Klee MS (2001) Quantitative comparison of performance of 
isothermal and temperature-programmed gas chromatography. J Chromatogr A 
933: 13-26. 

66. Castcllo G, Moretti P, Vezzani S (2009) Retention models for programmed gas 
chromatography. J Chromatogr A 1216: 1607-1623. 

67. Zhao C-X, Liang Y-Z, Fang H-Z, Li X-N (2005) Temperature-programmed 
retention indices for gas chromatography- mass spectroscopy analysis of plant 
essential oils. J Chromatogr A 1096: 76—85. 

68. Vogt T (2010) Phenylpropanoid biosynthesis. Mol Plant 3: 2-20. 

69. Giddings JG (1995) Sample dimensionality: a predictor of order-disorder in 
component peak distribution in multidimensional separation. J Chromatogr A 
703: 3-15. 

70. Schure MR (201 1) The dimensionality of chromatographic separations. 
J Chromatogr A 1218: 293-302. 

71. Oh C, Huang X, Regnier FE, Buck C, Zhang X (2008) Comprehensive two- 
dimensional gas chromatography/ time -of- flight mass spectrometry peak 
sorting algorithm. J Chromatogr A 1 179: 205-215. 

72. Wilson RB, Siegler WC, Hoggard JC, Fitz BD, Nadeau JS, et al. (2011) 
Achieving high peak capacity production for gas chromatography and 
comprehensive two-dimensional gas chromatography by minimizing off- 
column peak broadening. J Chromatogr A 1218: 3130—3139. 

73. Blumbcrg LM (2003) Comprehensive two-dimensional gas chromatography: 
metrics, potentials, limits. J Chromatogr A 985: 29-38. 

74. Dc Vos RC, Moco S, Lommcn A, Keurcntjes JJ, Brno RJ, et al. (2007) 
Untargetcd large-scale plant metabolomics using liquid chromatography 
coupled to mass spectrometry. Nat Protoc 2: 778-791. 

75. Xiao JF, Zhou B, Ressom HW (2012) Metabolite identification and 
quantitation in LC-MS/MS-bascd metabolomics. Trends Anal Chcm 32: 1- 
14. 

76. Grata E, BoccardJ, Guillarme D, Glauscr G, Garrupt P-A, et al. (2008) UPLC- 
TOF-MS for plant metabolomics: a sequential approach for wound marker 
analysis in Arabidopsis thaliana. J Chromatogr B 871: 261—270. 

77. Kuc J, RushJS (1985) Phytoalcxins. Arch Biochem Biophys 236: 455-472. 

78. Bateman KP, Castro-perez J, Wrona M, Shockcor JP, Yu K, et al. (2007) MS E 
with mass defect filtering for in vitro and in vivo metabolite identification. 
Rapid Commun Mass Spectrom 21: 1485-1496. 

79. Kind T, Fichn O (2007) Seven golden rules for heuristic filtering of molecular 
formulas obtained by accurate mass spectrometry. BMC Bioinfbrmatics 8: 1— 
20. 

80. Sumner LW, Amberg A, Barrett D, Beale MH, Beger R, et al. (2007) Proposed 
minimum reporting standards for chemical analysis. Metabolomics 3: 21 1-221. 



PLOS ONE | www.plosone.org 



17 



January 2014 | Volume 9 | Issue 1 | e87846 



Metabolomics of Ergosterol-lnduced Plant Responses 



81. Kanani HH, Klapa MI (2007) Data correction strategy for metabolomics 
analysis using gas chromatography-mass spectrometry. Mctab Eng 9: 39—51. 

82. Kanani H, Chrysanthopoulos PK, Klapa MI (2008) Standardizing GC-MS 
metabolomics. J Chromatogr B 871: 191-201. 

83. Ralston-Hooper K, Hopf A, Oh C, Zhang X, Adamec J, et al. (2008) 
Development of GCxGC/TOF-MS metabolomics for use in ecotoxicological 
studies with invertebrates. Aquat Toxicol 88: 48-52. 

84. Eriksson L, TryggJ, Wold S (2008) CV-ANOVA for significance testing of PLS 
and OPLS® models. J Chemometncs 22: 594-600. 

85. Sadeghi-bazargani H, Bangdiwala SI, Mohammad K (2011) Compared 
application of the new OPLS-DA statistical model versus partial least squares 
regression to manage large numbers of variables in an injury case-control study. 
Sci Res Essays 6: 4369-4377. 

86. Chong I-G, Jun C-H (2005) Performance of some variable selection methods 
when multicolinearity is present. Ghemometr Intell Lab 78: 103-112. 

87. Mehmood T, Liland KH, Snipcn L, Saebo S (2012) A review of variable 
selection methods in Partial Least Squares Regression. Chemometr Intell Lab 
118: 62-69. 

88. Trivedi KD, lies KR (2012) The application of SIMCA P+ in shotgun 
metabolomics analysis of ZIG HILIC-MS spectra of human urine - Experience 
with the Shimadzu IT-TOE and profiling solutions data extraction software. 
J Chromatogr Sep Tech 03: 1-5. 

89. Brown M, Wedge DC, Goodacre R, Kell DB, Baker PN, ct al. (2011) 
Automated workflows for accurate mass-based putative metabolite identifica- 
tion in LC/MS-dcrivcd metabolomic datascts. Bioinformatics 27: 1108—1112. 

90. Stoessl A, StothersJB, Ward EWB (1976) Sesquiterpenoid stress compounds of 
the Solanaceae. Phytochcmistry 15: 855-872. 

91. Chappcll J, Nable R (1987) Induction of sesquiterpenoid biosynthesis in 
tobacco cell suspension cultures by fungal elicitor. Plant Physiol 85: 469-473. 

92. Schnee C, Ko TG, GershenzonJ (2002) The maize gene terpene synthase 1 
encodes a sesquiterpene synthase catalyzing the formation of (E)-bcta- 
farnesene, (E)-nerolidol, and (E,E)-farnesol after herbivore damage. Plant 
Physiol 130: 2049-2060. 

93. McGarvey DJ, Groteau R (1995) Terpenoid metabolism. The Plant Cell 7: 
1015-1026. 

94. Nugroho LH, Pcltcnburg-Looman AMG, Verbernc MC, Verpoortc R (2002) 
Is accumulation of sesquiterpenoid phytoalexins induced in tobacco plants 
constitutively producing salicylic acid? Plant Sci 162: 989-993. 

95. Tholl D (2006) Terpene synthases and the regulation, diversity and biological 
roles of terpene metabolism. Curr Opin Plant Biol 9: 297-304. 

96. Chen F, Tholl D, Bohlmann J, Pichcrsky E (2011) The family of terpene 
synthases in plants: a mid-size family of genes for specialized metabolism that is 
highly diversified throughout the kingdom. Plant J 66: 212-229. 

97. Chappcll J, Von Lankcn C, Vogeli U, Bhatt P (1989) Sterol and 
sesquiterpenoid biosynthesis during a growth cycle of tobacco cell suspension 
cultures. Plant CeU Rep 8: 48-52. 

98. Zook M, Hohn T, Bonncn A, Tsuji J, Hammcrschmidt R (1996) 
Characterization of novel sesquiterpenoid biosynthesis in tobacco expressing 
a fungal sesquiterpene synthase. Plant Physiol 112: 311-318. 

99. Cane DE, Bowser TE (1999) Trichodiene synthase: mechanism-based 
inhibition of a sesquiterpene cyclase. Bioorg Med Chcm Lett 9: 1127-1 132. 

100. Oelofse D, Dubery IA (2006) Induction of defence responses in cultured 
tobacco cells by elieitors from Phytophthora nicotianae. Int J Biochcm Cell Biol 
28: 295-301. 

101. Hammcrschmidt R (1999) Phytoalexins: what have we learned after 60 years? 
Annu Rev Phytopathol 37: 285-306. 

102. Smith CJ (1996) Accumulation of phytoalexins: defence mechanism and 
stimulus response system. New Phytol 132: 1—45. 

103. Mert-Tiirk F (2002) Phytoalexins: defence or just a response to stress? J Cell 
Mol Biol 1: 1-6. 

104. Bari R, JonesJDG (2009) Role of plant hormones in plant defence responses. 
Plant Mol Biol 69: 473-488. 



105. Finkelstein RR, Rock CD (2002) The Arabidopsis book: Abscisic Acid 
biosynthesis and response. Rockville MD: American Society of Plant 
Biologists p. 1-52. 

106. Milborrow BV (2001) The pathway of biosynthesis of abscisic acid in vascular 
plants: a review of the present state of knowledge of ABA biosynthesis. J Exp 
Bot 52: 1145-1164. 

107. Seo M, Koshiba T (2002) Complex regulation of ABA biosynthesis in plants. 
Trends Plant Sci 7: 41-48. 

108. Anderson JP, Badruzsaufari E, Schcnk PM, Manners JM, Desmond OJ (2004) 
Antagonistic interaction between abscisic acid and jasmonate-ethylene 
signaling pathways modulates defence gene expression and disease resistance 
in Arabidopsis. The Plant Cell 16: 3460-3479. 

109. Mauch-Mani B, Maueh F (2005) The role of abscisic acid in plant-pathogen 
interactions. Curr Opin Plant Biol 8: 409-414. 

1 10. Zhang J, Jia W, Yang J, Ismail AM (2006) Role of ABA in integrating plant 
responses to drought and salt stresses. Field Crop Res 97: 1 1 1-1 19. 

111. Garcia-Mata C, Lamattina L (2007) Abscisic acid (ABA) inhibits light-induced 
stomatal opening through calcium- and nitric oxide-mediated signaling 
pathways. Nitric Oxide: Biol Chcm 17: 143—151. 

112. Toumi I, Mosehou PN, Paschalidis KA, Bouamama B, Salem-fnayou A Ben, et 
al. (2010) Abscisic acid signals reorientation of polyamine metabolism to 
orchestrate stress responses via the polyamine exodus pathway in grapevine. 
J Plant Physiol 167: 519-525. 

113. Song W, Ma X, Tan H, Zhou J (2011) Abscisic acid enhances resistance to 
Alternaria solani in tomato seedlings. Plant Physiol Biochem 49: 693-700. 

1 14. Mialoundama AS, Heintz D, Debayle D, Rahier A, Camara B, et al. (2009) 
Abscisic acid negatively regulates elicitor-induced synthesis of capsidiol in wild 
tobacco. Plant Physiol 150: 1556-1566. 

115. Jung HW, Tschaplinski TJ, Wang L, Glazebrook J, Grccnbcrg JT (2009) 
Priming in systemic plant immunity. Science 324: 89-91. 

1 16. Newman M-A, Von Roepenack-Lahaye E, Parr A, Daniels MJ, DowJM (2001) 
Induction of hydroxycinnamoyl-tyraminc conjugates in pepper by Xanthomo- 
nas campestris, a plant defence response activated by hrp gene-dependent and 
hrp gene-independent mechanisms. Mol Plant Microbe Int 14: 785—792. 

117. Kachroo A, Lapchyk L, Fukushige H, Hildebrand D, Klessig D (2003) 
Plastidial fatty acid signaling modulates salicylic acid-and jasmonic acid- 
mediated defence pathways in the Arabidopsis ssi2 mutant. The Plant Cell 15: 
2952-2965. 

118. Upchurch RG (2008) Fatty acid unsaturation, mobilization, and regulation in 
the response of plants to stress. Biotechnol Lett 30: 967-977. 

1 19. Kachroo A, Kachroo P (2009) Fatty acid-derived signals in plant defence. Annu 
Rev Phytopathol 47: 153-176. 

120. Wang K, Senthil-Kumar M, Ryu C-M, Kang L, Mysore KS (2012) 
Phytosterols play a key role in plant innate immunity against bacterial 
pathogens by regulating nutrient efflux into the apoplast. Plant Physiol 158: 
1789-1802. 

121. Madala NE, Molinaro A, Dubery IA (2011) Distinct carbohydrate and lipid- 
based molecular patterns within lipopolysaccharidcs from Burkholderia cepacia 
contribute to defence-associated differential gene expression in Arabidopsis 
thaliana. Innate Immun 18: 140-154. 

122. Gricbel T, ZeierJ (2010) A role for beta-sitosterol to stigmastcrol conversion in 
plant-pathogen interactions. Plant J 63: 254—268. 

123. Vogeli U, FreemanJW, GhappellJ (1990) Purification and characterization of 
an inducible sesquiterpene cyclase from clicitor-treated tobacco cell suspension 
cultures. Plant Physiol 93: 182-187. 

124. Keller H, Czernic P, Ponchet M, Ducrot PH, Back K, et al. (1998) 
Sesquiterpene cyclase is not a determining factor for elicitor- and pathogen- 
induced capsidiol accumulation in tobacco. Planta 205: 467—476. 

125. Granado J, Felix G, Boiler T (1995) Perception of fungal sterols in plants: 
subnanomolar concentrations of ergosterol elicit extracellular alkalinization in 
tomato cells. Plant Physiol 107: 485-490. 

126. Vatsa P, Chiltz A, Luini E, Vandelle E, Pugin A, et al. (201 1) Cytosolic calcium 
rises and related events in ergosterol-treated Nieotiana cells. Plant Physiol 
Biochem 49: 764-773. 



PLOS ONE | www.plosone.org 



18 



January 2014 | Volume 9 | Issue 1 | e87846 



